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I. INTRODUCTION 

Octupolar and tetrahedratic are synonymous adjectives when applied to soft matter ordering^ Loosely speaking, 
they are called upon whenever a third-rank tensor is needed to describe order in a molecular ensemble. Perhaps, 
Pel Hi was the first to consider a third-rank tensor to describe an unconventional liquid-crystal phase with the 
tetrahedron symmetry condensing directly from the isotropic phase0 While these early studies had concerned phases 
fully characterized by a third-rank tensor, a subsequent, comprehensive analysis considered the coexistence of polar, 
quadrupolar, and octupolar orders, elucidating the complex network of phase transitions thus made possible Q. Non¬ 
polar nematic phases, where only quadrupolar and octupolar orders can coexist, were studied in JO, where the reader 
is also referred to for a rather accurate and informative review of earlier contributions to theoryo Third-rank tensors 
do not only play a role in describing new condensed soft phases; they have also recently been employed in the active 
dynamics of self-propelled microorganisms. In a series of papers extending the ideas originally presented in 

[nl, a third-rank tensor is invoked to represent the asymmetric shape of a living cell when spontaneous deformation 
and drift velocity are intimately interconnected. 

Here our attention is mainly directed to describing octupolar order through a third-rank, fully symmetric and 
traceless tensor. It is known that in three space dimensions that are 7 independent such tensors they become 
2 in two space dimensions. To diagonalize a third-rank tensor is not a univocally defined task, as no analogue of 
the Spectral Theorem exists for symmetric tensors with rank higher than 2. Likewise, a commonly accepted notion 
of generalized eigenvalues and eigenvectors has not been established yet for these tensors, which would univocally 
identify the former as scalar order parameters!! We shall find a viable substitute for such a classical representation of 
the scalar order parameters in the maxima (and conjugated minima) of the octupolar probability density that shall 
be introduced in Sec. im 

In this paper, only the two-dimensional case will be considered: it is far easier than the three-dimensional case (which 
will be studied elsewhere El), but not unrealistic, as suggested by recent experiments with nails on a vibrating table 
Ellil- Compared to either tetrahedratic phases or self-propelling microorganisms, this is perhaps a more mundane 
manifestation of the need for an octupolar order descriptor. Section IIIII is a brief interlude with the purpose of 
reinterpreting the traditional descriptions of both dipolar and quadrupolar orders in terms of maxima (and conjugated 
minima) of the appropriate probability density. In Sec. IIVI which is the heart of the paper, the octupolar order in 
two space dimensions is formally characterized in terms of the generalized eigenvalues and eigenvectors of a third- 
rank tensor, which is also given an equivalent diagonal representation, whose validity is however restricted to the 
two-dimensional case. In Sec. El both strategy and conclusion of this paper are recapitulated and their foreseeable 
extensions to the three-dimensional case are briefly anticipated. 


II. PROBABILITY DENSITY MULTIPOLES 


Consider, for definiteness, a system of rigid nails in a plane like those studied experimentally on a vibrating table 
TtI . Each nail is characterized by a unit vector p oriented from its head to its end (see Fig.[T]). Assume that nails 



FIG. 1. (Color online) Cartoon illustrating a nail and the polar unit vector p represented as in (I13II relative the Cartesian 
frame (£ 1 , 62 ). 

are distributed in the plane with a certain probability law which may induce order in their ensemble. The probability 


^ Despite their equivalence, we shall use throughout the former instead of the latter. 

^ A wealth of unconventional nematic phases allowed by symmetry are described in Q. Pel’s analysis was revisited in Q. 

® A full macroscopic theory is also proposed in Q, which encompasses statics and hydrodynamics of these phases, which are referred to 
as D2d phases. In particular, a macroscopic dynamic theory for the octupolar T^i phase is presented in Sect. 2.2 of Q, building on an 
earlier work Q. 

This issue has already been addressed in a previous study [T3l which proposed an approach alternative to the one followed here. Likewise, 
other computational definitions of scalar order parameters for both tetrahedral and cubatic symmetries can also be found in [Ulil- 
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distribution density, which is defined on the unit circle will be denoted by 
normalization condition 


/SI 


g{p)ds{p) = 1 , 


■; it is subject to the 


( 1 ) 


where s{p) is the arch-length on Adapting the formalism of [T^ to represent Buckingham’s formula [T^ for g in 
terms of Cartesian tensors, we write 


9ip) = ^ 1 + 


f=i 


where is the j-th rank tensor defined by 


:= p ( 




( 2 ) 


( 3 ) 


j times 


C) denotes tensor product, and is a jth-rank tensor soon to be related to the jth moment of pH In ([2), the 
inner product of two tensors, and of equal rank j, corresponds to the following contraction of 

components in any Cartesian frame ( 61 , 62 ) of the plane: 


A« -BW := 


( 4 ) 


Here and in what follows we understand that repeated indices are summed over their whole range. 

Since the products in ([3]) are completely symmetric under the exchange of any pair of components for all i. by (|3|) 
so are also required to be all P^-^^ to dispose of redundant components. For the function in ([3]) to obey ([T]), all P^-') 
are further required to be traceless in any pair of components, as isotropy of the plane demands that the tensor 


■= 


2it 


p^^ds{p) 


( 5 ) 


either vanishes, if j is odd, or its Cartesian components can be written as symmetrized products of Kronecker’s 5’s, 
if j is even. Letting the brackets TTT' denote the irreducible, completely symmetric and traceless part of any tensor 
they surmount, the above properties of each are embodied by the equation 


pU) = p(f) V j 


We shall denote by (...)„ the ensemble average relative to g: 


(...) := f {...)g{p)ds{p), 


( 6 ) 


( 7 ) 


so that the average (.. .)q in ([5]) corresponds to the average relative to the isotropic density function go = 
It is a direct consequence of ([5]) and (HI) that 

p'o’j ^ o P^'^^ V j ^ I, 

where o denotes tensor multiplication^ While it is almost immediate to prove from (|S]) that 

1 


_ 1 


it requires some tedious labour to show that 


P 


iP), = 2 ^^'^ 


= _p0) V 7 > 2, 

27 7 - > 


( 8 ) 


(9) 


( 10 ) 


^ In general, q could alternatively be represented as an expansion in symmetry-adapted Wigner rotation matrices |[2Q|| . but both here and 
in a forthcoming paper [Tsil we are interested in the equivalent Cartesian tensor representation of q. 

® If the Cartesian components of tensors A and B, of rank 2j and j respectively, then C = A o B is a 

tensor of rank j and its components are given by 
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so that we can rewrite © as 


eip) 



( 11 ) 


In (|lll) . g is expressed as the sum of density multipoles, each associated with a corresponding order tensor0 We are 
especially interested in the first three order tensors, {p)g, , and , featuring in (ITTI) . which 

represent three independent descriptors of order; we call them dipole, quadrupole, and octupole, respectively. The first 
and the last are measures of polarity, the latter becoming relevant when the former is bound not to dominate, as 
suggested by the tendency of shacked nails shown in [l^, [l^ to be on average combined in antiparallel pairs. For a 
pair of nails, it was indeed shown in (2lj | that in the antiparallel configuration the excluded volume is smaller than in 
the parallel configuration!! 

The strategy pursued here is to identify the scalar order parameters of the order tensors in (fTTl) with the maxima 
(and conjugated minima) of the density multipoles. From now on, to avoid clutter, we shall drop the subscript g from 
the averages (...)^, as g will be the only probability density we shall consider in the following. 


III. DIPOLAR AND QUADRUPOLAR ORDERS 

Although this paper is concerned with octupolar order in two space dimensions, I find it useful to indulge in 
rephrasing both dipolar and quadrupolar orders in terms of maxima (and conjugated minima) of the appropriate 
probability density multipole. 


A. Dipole 

The average dipole (p) is a vector in the plane, which can be represented as 

(p) = Aid, (12) 


where d is a unit vector and Ai, which can be taken as positive, is the dipolar scalar order parameter. Letting 
p ■ d = cosd, it follows from (fT^ that Ai = (cosd) and so 0 ^ Ai ^ 1. Alternatively, representing p in a Cartesian 
frame ( 61 , 62 ), as in Fig. [U 



p = cos 62 + sin ipei. 

(13) 

and setting 

d := costpo 62 + sinipo 61 , 

(14) 

we see that 

Ai sin ipQ = (sin ip) and Ai cos po = (cos p). 

(15) 

Then the function 

Pi{p) ■= -{p) -P 

TT 

(16) 

expresses the dipolar density in 

pip. Once normalized to its maximum, it simply becomes 



Pi (p) = d • p = cos{p - pq) 

(17) 


and its polar plotj! shown in Fig. [U is a circle passing through the origin. It should be noted that for pi(p), as for 
any odd function on the polar plot is actually drawn twice as p ranges over S^: in one wrapping, pi is positive, 
whereas it is negative in the other. Thus, the polar plot of pi is just the same as the polar plot of its positive part 
/x)" := max{0,pi} (or its negative part pj" min{0,pi}, for that matter). Clearly, as shown in Fig. |2l the dipole d 
corresponds to the direction of p with maximum dipolar density, Ai/tt, whereas the direction — d is corresponds to 
the direction of minimum dipolar density, —Ai/tt. 


^ Which represents density moments of the appropriate rank. 

® For cylindrically symmetric, convex bodies, possibly tapered along the symmetry axis, it has recently been shown that the antiparallel 
configuration actually minimizes the excluded volume [ 2 ^ , a result that had been proved earlier for cones 1231 . How general this property 
could be for cylindrically, non-convex shapes is still a matter of debate. 

® The polar plot a scalar function p(p) is the curve described on the plane by the vector p(p)p as p ranges over . 
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FIG. 2. (Color online) The dipole d superimposed to the polar plot of the normalized dipolar density pi in (1171) . which illustrates 
pictorially the dipolar probability density. Here ipo = f. The vector —d, which represents the least probable orientation of p, 
is dashed. 


B. Quadrupole 


Similarly, the quadrupolar order tensor 
represented as 


p®p 


is a second-rank tensor which by the Spectral Theorem can be 


p®p) = \2e®e, 


(18) 


I- 1 


where e is the unit eigenvector of i^p® p j associated with the positive eigenvalue A 2 . Letting now p ■ e = cos'd, it 
follows from (fTSl) that 

(cos^d) = A2, (19) 

and so 0 ^ A 2 ^ 1. SettinsF^ 

e = sin (po ei-b cos (po ^ 2 , (20) 

from (fT^ we also obtain that 

A 2 sin 2 (po = (sin2(p) and A 2 cos 2 ipo = (cos2(p). (21) 

The mapping defined on by 

2 


P2{p) := - (^p®p'^ ■ p® p 


( 22 ) 


expresses the quadrupolar density in CJ). In complete analogy with pi in mi), we write the normalized quadrupolar 
density as 


P2(j 5) = 2 e 0 e • p (g) p = cos(2(p — 2(po). (23) 

The polar plot of the positive part ^ = max{0,p2} of p 2 is depicted in Fig. [3] along with the unit vector e in (1301) 
and its orthonormal companion e±. The quadrupolar density p 2 has equal maxima, A2/'7r, along e and — e and equal 
minima, —A2/7r, along and —ej_. 

The rationale behind plotting only the positive part pj, which is our choice here, is the desire of characterizing the 
probability density multipoles (at least in two space dimensions) in terms of their maxima (and conjugated minima) 
and the directions in where they are attained. As will appear clearer in our development below, this is in tune 
with the notion of generalized eigenvalues and eigenvectors that shall be employed to describe the octupolar order (as 
well as possibly also higher-rank orders). 


It is perhaps worth mentioning that the angle (po designating e through 112011 and featuring in dn]) need not be the same as the angle 
denoted in the same way but featuring in Gl, as there is no guarantee that d and e should either coincide or be somehow related. This 
slight abuse of notation is not likely to confuse the reader, if one heeds that only d and e are physically relevant and not the angles that 
designate them in the plane. 
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FIG. 3. (Color online) The polar plot of pj = max{0, P 2 } with p 2 as in (I23II . Here po = f • Both unit vectors e± and —e± are 
dashed. 


IV. OCTUPOLAR ORDER 


Representing the octupolar order tensor 
form 


I-1 

p'Sip'Sip 


is not as simple as writing the analogue of (fT^ in the diagonal 


pSpSp) = X'ia 


' a S a, ■ 


(24) 


as we lack a generally accepted notion of eigenvalues and eigenvectors for tensors of rank higher than 2 , and, more 
importantly, for such tensors we lack the analogue of the Spectral Theorem. We shall see below that (l24)) can 
indeed be justifiedlli] but this requires some labor and resort to an appropriate notion of generalized eigenvalues and 
eigenvectors. Here we start from the Cartesian representation 

(^pSpSp'j = athkeiS ehS ek, (25) 

where 

a^hk ■■= (piPhPk - ^ {PiShk +Ph5ik +PkSih)^ , (26) 

having denoted by pi the components of p in the Cartesian frame (ei, 62 ). It readily follows from (l26l) that 


Oiii + 0122 = 0, 

0211 + 0222 = 0 , 

— Olll = 0122 = O22I = O212, 

— 0222 = O211 = O112 = O121, 


(27) 


which show that only two components aihk are indeed independent; we select oi := om and 02 := 0222 to represent 
all of them. By using again m, we obtain from (|26l) that 

oiii = — - (sin3(p) and 0222 = ^ (cos3(p). (28) 

In accord with our earlier treatment of both the dipolar and quadrupolar densities in (HU, we effectively represent 
(^p S p S p^ through the maxima (and conjugated, opposite minima) of the function pa defined over by 


P3{p) ■■= ^(^pSpSp'j 


■pSpSp, 


(29) 


which designates the octupolar density. Looking for the constrained extrema of ps over amounts to solving the 
problem 


^ihk^h^k ^ 


(30a) 
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In a way that actually makes the two-dimensional case exceptional. 
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for a; G ad A G M, where A is the Lagrange multiplier associated with the constraint 

xf + X 2 = 1- (30b) 

The solutions (A, a:) of problem (l30l) coincide with the generalized eigenvalues and eigenvectors of (^p iS> p p'^ 

according to a definition introduced for tensors (not necessarily symmetric) of rank higher than 2 in finite-dimensional 
spaces of any dimension. This concept has been proposed and made precise in (23 - [^ . Different notions of generalized 
eigenvalues and eigenvectors have been introduced in the literature; for the one chosen here, a theorem was recently 
proved in on the cardinality of the eigenvalues, which to my knowledge is not available for other notions. According 
to this theorem, in the case at hand there are at most three distinct complex eigenvalues (defined to within a sign), 

of which at least one is real. We shall show below that the solutions of (1301) are rather easy to find: they turn out to 

be all real and equal to one another, to within a sign. 

By use of (Eli), we write (I30al) in the explicit form 

ai{xi — X 2 ) — 2a2XiX2 = Axi, (31a) 

02 ( 0:2 — xf) — 2aiXiX2 = Ax 2 , (31b) 

which make it evident that to each solution (A, a;) there corresponds the solution (—A,—a:). To solve (1311) subject to 
(I30bl) we oarameterize the latter bv lettine. as in dSD, 

X 2 = COS If and xi = sin(/j. (32) 


The parity symmetry (A, a:) >->• {—X,—x) thus translates into {—X,(p + tt), so that each solution (p G [0,7r] 

generates another solution in [tt, 27r] by a 7r-shift. Inserting (15^ into (I31|) . we arrive at 

—oi cos 2(^ — 02 sin 2(^ = A sin (33a) 

02 cos 2(^ — oi sin 2(^ = A cos (33b) 


We distinguish two cases, oi = 0 and oi ^ 0. In the former case, eliminating A from (1551) we obtain 

tan 2p + tan p = 0, 

which in [0, tt] has solutions pi = 0, P 2 = f , and ps = |7r. Correspondingly, A is delivered by A*^^) = 02 , A(2) 
and A(3) = 02 . If oi ^ 0, we set a := ^ and eliminating again A we obtain 

1 

•a =-) 

tan 6p 


(34) 


- 02 , 


(35) 


which has three solutions in [0, tt], denoted pi, p 2 , and ps, to which there correspond three values of A, A*^®^ = aiA{pi), 
delivered by the function 


A((p) := 


1 

sin3(/? 


(36) 


The strategy to solve (1551) for oi 7 ^ 0 is illustrated graphically in Fig. |4l For a generic a, the pi’s are obtained by 
intersecting the line y = a with the graph of y = — l/tan 3 (/ 5 . Correspondingly, the A^*^ (scaled to oi) are read off 
from the graph of y = A(y)). Since A((^ ± ^) = —A(yj), and any two adjacent roots pi differ by the corresponding 
A*^*^ alternate in sign and for all of them |A^ = loilvT -f . The solution shown in Fig. |4]is for a > 0 and delivers 

A^^) < 0 if also oi > 0. For a < 0 and oi > 0, A^^l would also change its sign. 

It is worth noting that all extrema of pa on are either maxima or minima (conjugated by parity to the former): 

no saddle thus corresponds to a generalized eigenpair of (^p ® p ® p^. The eigenvectors ai = x^^ eu of (^p®p®p^ 

are then 6, counting all solutions of (1311) conjugated by parity. They will be conventionally represented by solid 
vectors if corresponding to positive eigenvalues and by dashed vectors if corresponding to negative eigenvalues. The 
solutions shown in Fig. |5] are reproduced in Fig. [5] with this convention. 

Rescaling the octupolar density p?. in (|2^ to its maximum, we reduce it to a simple function of p, pa = cos(3(p—3po)j 
where po is any solution of (|55|) Fd Compared to the polar plots of pi and pj in (fT71) and ((55|) . which exhibit one and 


Again, I am guilty here of a slight abuse of notation, as (/jq IIS'S already been used above to represent d and e, and the corresponding 
normalized densities pi and p 2 . No relation should be expected between these angles, though they are occasionally denoted in the same 
way. 














FIG. 4. (Color online) Graphical illustration of the solutions of (1331) for m 7 ^ 0. The </5i’s, which are delivered by the 
intersections (hollow circles) oi y = a (solid thin line) with y = —Ij tan3v5 (solid thick lines), differ by ^ from one another. 
The corresponding which are read off from the graph of A(ip) = —1/sin 3^3 (full circles on the thick dashed lines), alternate 
in sign. Here 0=1 and |A^‘^| = ||cn|. 
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FIG. 5. (Color online) The 6 unit vectors ai, conjugated by parity, that correspond to the solutions of (1311) shown in Fig. |4] 
The dashed vectors ai, — 02 , and 03 , which correspond to the eigenvectors of the tensor in (1251) with negative generalized 
eigenvalues, represent the orientations of p with the least octupolar probability density ps. The superimposed polar plot 
represents the positive part of the normalized octupolar density ps. 


two lobes respectively, the polar plot of ^ := max{ 0 ,p 3 }, the positive part of pa, displays one more, as expected; in 
Fig. El it is superimposed on the 6 dials that represent the eigenvectors a^. 

Here we have chosen to characterize the third-rank tensor in terms of its generalized eigenvalues and 

eigenvectors corresponding to maxima and minima of pa. However, in the two-dimensional setting, equation p4|) is 
also proved valid letting a be any generalized eigenvector of and Aa = 4A, where A is the eigenvalue 

associated with the selected a. To see this, it suffices to represent a in (l24l) as 

a = cos (p 62 -I- sin (p ei, (37) 

which leads us to identify the two independent components of as 

oi = Aa sin tp (sin^ p — |) , (38a) 

02 = Aa cos p (cos^ p — |) . (38b) 

Solving these equations for A 3 and p, we readily obtain that A 3 = 4aiA(p), where A(p) is as in (j55)) . and p is a root 
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of (1551) . We thus conclude that in the two-dimensional case the generalized eigenpairs (A,a) of afford 

a 6-fold degenerate representation of this third-rank tensor through (IMl) . in complete analogy to what the Spectral 
Theorem does for a second-rank tensor[^ It will be shown in [l5l| that this is indeed a rather exceptional circumstance. 


V. CONCLUSIONS 


The strategy inspiring our quest for the representation of the octupolar order in any space dimension can be easily 
summarized by saying that no diagonalization need in general be attempted for a third-rank fully symmetric and 
traceless tensor (an octupolar tensor, for short). This latter should rather be characterized in terms of the maxima 
(and conjugated, opposite minima) of the octupolar probability density. In algebraic terms, this amounts to compute 
the relevant generalized eigenvalues and eigenvectors according to a definition for which a theorem concerning their 
cardinality has recently been proved 

In two space dimensions, an octupolar tensor is simply described by 2 scalar parameters, and so its generalized 
eigenvalues and eigenvectors in the plane must be so constrained as to be described by 2 parameters only. It was 
indeed shown here that in the two-dimensional case all generalized eigenvalues of an octupolar tensor are equal (to 
within a sign) and that its three inequivalent eigenvectors are described by a single rotation angle. It will be shown in 
M how the 7 independent parameters that describe an octupolar tensor in three space dimensions concur to either 
the three or four inequivalent generalized eigenvalues and eigenvectors which are associated in the generic case with 
the maxima of the octupolar probability density and the directions along which they are attained. 

A critical voice could be raised against the mathematical machinery devised here to describe the multipoles of the 
probability density g in (ED. Comparing pi, p 2 , and ps above 0 one could easily argue that the multipoles in (ITTI) are 
nothing but the Fourier components of p in . Such a critique could indeed be grounded if the angle ipo were indeed 
the same in all the formulae where it appears, which would the case if the vectors d, e, and a were somehow locked, 
that is, rigidly related to one another, which they are not. In general, there is no reason (apart from convenience and 
laziness) why one should think that dipolar, quadrupolar, and octupolar order tensors have correlated eigenvectors. 

In two space dimensions, a diagonalized form for an octupolar tensor was established in (1241) . which unfortunately 
has no analogue in three space dimensions 1^. It is thus conceivable that the method proposed here to identify the 
scalar order parameters of an octupolar tensor could be extended to any space dimension, whereas the diagonal form 
in (IMl) is only accidentally valid in two space dimensions. 
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